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BACKGROUND OF THE INVENTION 

rOOOn 1. Field of tho Inv o ntion SYSTEM FOR DUST DETECTION AND 
DUST IMAGE ENHANCEMENT 

[00011 BACKGROUND OF THE INVENTION 
100021 Field of the Invention 

[0002] This [0003| This invention deals generally with digital enhancement of 
e nvironm e ntal s at e llit e imag e ry of airborn e dust layers ov e r both oe e on and land 

backgrounds, meteorological imagery of dust clouds over ocean backgrounds. 

100031 2. D o Gcription 100041 Description of the Related Prior Art 

[000 4 ] Th e [00051 The detection and tracking of airborne dust from satellite has been 
efa long-standing interest to both gf the Navy and academia alike, and there is no shortage 
of papers treating the subject (e.g., Shenk. W. E., and R. J. Curran. 1974: The detection of 
dust storais over land and water with satellite visible and infrared measurements, Mon. Weather 
Rev., 102. 830-837. Carlson, T. N., 1978: Atmospheric turbidity in Saharan dust outbreaks as 
determined by analyses of satellite brightness data. Moii. IVeatlwr Rev.. 107, 322-335. Legrand, 
M., M. Desbois, and K. Vovor. 1987: Satellite detection of Saharan dust: optimized imaging 
during nighttime. J. Climate. 1. 256-264. Jankowiak, I., and D. Tanre, 1992: Satellite 
climatology of Saharan dust outbreaks: method and preliminary results. J. Climate. 5. 646-656.) 



^[00061 



These earlier works f 



contrast and spectral contrast signatures to identify and retrieve properties of dust against 
land and ocean backgrounds. 



rOOOSI A fOOOTI 



A recent study (R. Cantu; ^^The 'The Role of Weather in Major 



Naval Aviation MishapG". Mishaps'\ MS Thesis, 106 pp., NPS Monterey) fiadshas 
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indicated that poor visibility has factored into approximately 54% of Navy aircraft 
"Class - A" "Class-A" mishaps, for an estimated annual loss of S51 million. Of these 
mishaps it was reported that 56% were preventable had better forecasting and/or 
obserx'ational tools been available. The current Dust Enhancement Product represents a 
paradigm shift in the way dust is observed from the space platform, and is directly 
applicable to mitigating this multi-million dollar liability. 

[00061 As noted in the roferenc e listing above, oovoral altomativo l 00081 Alternative 
methods exist for enhancing dust over land and water. Single-channel 
radiom e t e r imagery with scaling thresholds chosen to enhance a small dynamic range of 
temperature or reflectance values has been applied to visible and infrared geostationary 
imagery. The resultant enhancements, usually only of marginal visual quality, tseis 
highly dependent on the tliresholds chosen and hence not suitable for operations where 
optimal values for these tliresholds may change dynamically over space and time. 
Combinations of single-channel visible and infrared data is a marked improvement over 
the scaling described above, but still suffers cloud/dust ambiguity for lack of multi- 
channel visible data (i.e., the basis for the over-water Dust Enhancement Product, which 
requires blue and short wave infrared channels). The Empirical Orthogonal Function 
(EOF, or sometimes referred to as "principle componcnt") "principle component") method 
has been shown to do a good job of enhancing dust over land and water. However, results 
produced by this method again are highly dependent on the variability of the scene and 
hence will provide an inconsistent enhancement depending on background and dimension 
of the region. Spatial contrast reduction, which takes advantage of the blurring 
e^e taffect of dust over otherwise sharply- sharply defined lenestrial features, has been 
applied successfiilly for dust detection over land, but eften naturally fails over laminar 
ocean backgrounds and is limited over land in regions wherethe terrain is less variable 
(e.g., desert plains, where dust incidentally is most common). 

[0007] Non e of th e pr e vious m e thods enlist f00091 None of the above-mentioned 
methods enlists multi-channel visible data to take advantage of the inherent spectral 
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variability of dust at these wav e l e ngths, du e primarily to the unavailability of such data at 
the tim e s of those wTitings. wavelengths. The unified land/ocean Dust Enhancement 
Product is a Heve lunique solution to the daytime dust detection problem, fe 
comparisons A comparison between the visible/infrared combination technique (e.g.. as is 
applie d(applied here to Meteosat 5 data) and the current algorithm (applied to Terra- 
MODIS data) for a space and time co-located dust storm event in Southwest Asia has 
r e v e al e d been conducted. The superior detection of dust by the current method over both 
land and wat e r. As such, the n e w tochniqu e is water as demonstrated by these data are 
thought to be of immediate relevance to Navy Meteorology/Oceanography (METOC) 
operations in several regions of the worid that experience significant dust outbreaks 
(including Southwest Asia, Africa, the Mediterranean, and the West Pacific). 

SUMMARY OF THE I>rVENTION 100101 SUMMARY OF THE INVENTION 
[0008] The fOOUIThe dust enhancement system of the invention is dcsigned to 
provide an indication of the horizontal distribution of airborne dust lay e rs over both 
land and_water surfaces for daytime scenes. The regions of dust are enhanced as shades 
of pink, depending on opacity, against darker non-dust backgrounds (e.g., 

blu e wat e r, and cyan clouds). 

[00091 Tho [0012] The primary advantage of the Dust Enhancement Product is its 
ability to call immediate attention to areas where s ignificant concontrations o ft here 
airborne dust resid e , particulariv residcs. often in regions where even the true color 
imagery experiencesgfeat difficulty in discerning. The enhancement highlights 
mesoscale circulations in optically thin dust ov e r water , and is also useful for identifying 
volcanic plumes (e.g., Mt. Etna Jul. 26, 2001 e ruption imago 7/26/01 eruption — image 
attached). The latter capability is of particular relevance to both military and civilian 
aviation. The algorithm is applicable to a wide range of ocean color instruments in the 
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ocean-only mode (i.e., true color over land surfaces) or in the ocean/land mode provided 
infrared information is available either on the same platform or co-registered with 
temporally matched imagery from an independent observing system. 

[0010] The now dust enhanc e m e nt t e chniqu e is in fact a combination of two ind e p e nd e nt 
enhancement algorithms on e applicabl e to water backgrounds and the other to land 
backgrounds. It i s bas e d physically on th e pr e mise that silicat e a e rosols (e.g., mineral 
dust common to d e s e rt sandstorms) pr e f e r e ntially absorb blu e light compared to groon 
and r e d light, wher e as oloud particl e s (liquid or ice) are more or less non absorbing 
across this r e gion and refl e ct all thre e compon e nts strongly. Therefore, a sp e ctral 
diff e r e nce betw ee n r e d and blue light will b e much larg e r for dust than for clouds. In this 
way, data from a radiometer with sp e ctral r e solution capable of m e asuring indep e ndently 
th e r e d. gre e n, and blu e components of r e fl e ct e d solar energy can be combined to further 
enhanc e and decouple the dust signal from clouds. Over water, the spectral diff e renc e is 
normalized by the magnitud e of the reflectanc e such that additional s e nsitivity to thin 
(low r e flcotanco) dust is achiev e d. The approach r e quir e s correction of the multi channel 
data to r e move scatt e ring contamination from the molecular atmosph e re prior to 
computing the enhancem e nt. 

[0011] B e cause the land backgrotmd poss e ss e s similar spootral prop e rties to the airborn e 
dust at visiblo/shortwa^^e infrared wav e l e ngths, the simple oc e an component describ e d 
above is insuffici e nt for identifying dust ov e r land. Th e now toohniquo includ e s th e rmal 
contrast and spectral fransmittanc e prop e rty diff e r e nc e s, which requir e chann e ls in th e 
infrar e d part of tho sp e ctrum. With its 36 narrowband chann e ls spaiming the optical 
portion of the e lootromagnetic sp e ctrum fro 0. 4 to 1 4 . 4 micromet e rs, MODIS is w e ll 
e quipped to exploit those tochniquos from a single platform (although tho Dust 
Enhancoment Product may b e achi e v e d by combining r e gister e d data from a coUoction of 
ind e pend e nt platforms). 
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[0012] The three factors contributing to th e ov e rland compon e nt of th e Dust 
Enhanc e ment Productor e outlin e d as follows: During the day, cool e r e missions from an 
el e vat e d dust lay e r contrast against th e warm e missions from a h e ated surface. Th e some 
t e mp e rature signature holds tru e for clouds, but th e sam e sp e ctral d e coupling applied to 
th e over ocean e nhancem e nt appli e s also to clouds ov e r land. Combining th e t e mperature 
and shortwav e differenc e s provide s a m e ans to d e t e cting th e optically thick portion of the 
dust plume. Optically thin r e gions of dust will not produc e a strong infrar e d t e mperatur e 
d e pr e ssion. — Howev e r, — sp e ctral — differences — m — th e rmal — (11.0 12.0 — microm e t e rs) 
transmission for dust ar e very useful in determining ar e as of thin dust. The signature is 
opposit e in sign to that of thin cirrus. Combining these multisp e ctral f e atures together 
yi e lds a tractabl e algorithm for th e e nhancem e nt of dust over land. Although th e land and 
wat e r e nhanc e ment algorithms diff e r significantly, the land e nhanc e ment has been 
formulated such that dust over land is the sam e hu e (shades of pinlc) as dust ov e r wat e r 
with minimal discontinuity in saturation across coastal boundari e s. 

[00131 The [00131 The techniques of the invention demonstrate that the current 
onhancomont. techniQue. as applied here to the Sea-viewing Wide Field-of-view Sensor 
(SeaWiFS) radiom e ter (ov e r wat e r only) and MODIS (water and land). radiometer. has 
significant potential to benefit the warfighter through its ability to provide through high 
spatial resolution satellite imagery a detailed depiction of atmospheric dust outbreaks 
with enhancements for optically thin regunes, as well as other meteorological/land- 
surface features of interest (e.g., clouds, fires, snow cover). Based on these findings, 
NRL and Fleet Numerical Meteorology and Oceanography Center (FNMOC) have 
collaborated to host a subset of these products upon the FNMOC Secure Internet 
(SIPRNET) beta-webpagc for the purpose of i) making this new technology available to 
Navy assets, and ii) receiving feedback usefiil for improving upon its current utility. 

[00141 BRIEF DESCRIPTION OF THE DRAWINGS 

rOOHl FIG. 1 Tho lOOlSI Fig. 1 illustrates a SeaWiFS sensor (left) flying aboard the 
SeaStar platform (right) (right). 
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i y[00161 Fig. 2 illustrates raw SeaWiFS data from a Navy Regional 
Center (NRC) receiving station (left) and the noise/line corrected equivalent (right) 
(right). 

[00161 FIG. 3 SoaWiFS rOOni Fig. 3 illustrates a SeaWiFS High Resolution 

Picture Transmission (HRPT) tnie color image of the Eastern Atlantic Ocean, Iberian 
Peninsula, and Northwestern Africa, captured by NRC in Rota, Spain (F e b. 13, 
2004^ (02/13/01. 1255Z). A large dust plume wrapping into a baroclinic system is 
indicat e d indicated. 

[00171 FIG. 1 ScaWiFS fOOlS] Fig. 4 illustrates a SeaWiFS HRPT true color image 
of the Persian Gulf area, captured by NRC in Bahrain. Airborne dust is indicated. 
Phyloplankton is responsible for the green pigmentation in the littoral regions regions. 

rOOlSIFIG.S SoaWiFS r00191 Fig. 5 illustrates a SeaWiFS HRPT true color image 
of the Korean Peninsula, captured by NRC in Yokosuka Japan. Yellow arrows indicate 
areas of biomass s mok e /pollution smoke/pollution. 

[00191 FIG. 6 - - Truo f 00201 Fig. 6 illustrates a true color (top) and vegetation 
enhancement (bottom). Bright red indicates regions with abundant plant life. Note 
seasonal changes in Fall vegetation compared against Summer (FIG. 5) ( Fig. 5). 

[00201 FIG. 7 Coooptual illustration o fT 00211 Fig.7 illustrates preferential 
absorption of blue light by mineral dust, resulting in yellow reflected light. Clouds 
scatter more uniformly over the visible wavelengths, resulting in white t^fe Might. 

[00211 FIG. 8 Exainplcs [00221 Fig. 8 illustrates examples of thin dust enhancements over 
the waters of the Persian Gulf. Optically thin dust plumes appear as darker tones of red. 
Land regions in the dust enhancement have been masfee^- masked. 

[00221 FIG. 9 -- Tru e |00231 Fig. 9 illustrates true color imagery of Typhoon Pabuk 
making landfall on Japan Japan. 
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[00231 FIG. 10 Th e [0024| Fig. 10 illustrates an enliancement captur e s a captured 
volcanic ash plume from the 2001 emption of Mt. Etna. Land regions in the dust 
enhancement have been masked eat-out. 

[00211 FIG. 1 1 MODIS f00251 Fie. 11 illustrates a MODIS dust enhancement (0.5 
km resolution) over the Northern Arabian Sea, true color over land. Dust over the ocean 
and appear as shades of pink. Dry lakebeds over land appear as patches of pink. 
Accompanying text from Navy MetOc post deployment report corresponding to the 
storm observed in this image. 



[00261 DETAILED DESCRIPTION OF THE INVENTION 

[00251 The |0027| The Satellite Meteorological Applications Section at the Naval 
Research Laboratory (NRL) in Monterey has developed tme color, vegetation 
enhancement, and airborne sand/dust enhancement processing capabilities for the Sea- 
viewing Wide Field of View Sensor (SeaWiFS) instrument. These products are made in 
near real time (typically 60 to 90 minute latency from time of ingest at the receiving 
stations) using telemetries captured at Na\7 Regional Centers (NRCs) located in Spain 
(Rota), Bahrain, and Japan (Yokosuka), and populate a password-protected website for 
several predefined sectors within the general coverage regions of these centers. 
Presented here is an overview of the SeaWiFS sensor and its capabilities/limitations, a 
summary of the procedure, and the current processing architecture. Examples drawn 
from the archived database (automated processing of these data commenced Attfr August 
8, 2001) are presented to illustrate their quality and potential operational utility. 
[00261 Platform/InGtrum e nt Specifications utility. 
[00281 Platform/Instrument Specifications 

[00271 Tho l00291 The SeaWiFS instrument (Firestone and Hooker, 1996), launched 
aboard an extended Pegasus rocket on AagrAugust 1, 1997, is an 8-channel whiskbroom 
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scanning radiometer featuring 8 narrowband (. about. 20 (~20 nanometers (nm) wide) 
spectral channels centered at 412, 443, 490, 510, 555, 670, 765, and 865 nm wavelengths. 
Carried aboard the SeaStar spacecraft (see FlGrFig, 1), SeaWiFS follows a sun- 
synchronous (9 8 .3.degr e e. f98.3° inclination, 705 kilometer (km) altitude above mean sea 
level (AMSL)) orbit with a local noon descending equatorial crossing time. It has a 
local-area-coverage (LAC) swath width of 2801 km (corresponding to a maximum scan 
angle of 5 8 .3.d e gr ee .) 58.3°) and a nominal spatial resolution of 1.1 km at nadir. Global 
Area Coverage (GAC) is also available at 4 km degraded resolution and a reduced 1502 
km swath (data limited to 4 5.degr ee . 45° scan angle). The 10-bit digitized data are 
transmitted dowTi to Earth at an L-band frequency of 1.70256 Gigahertz (GHz). The 
instrument scans across 1285 picture elements at a rate of 6 Hz. Typical High Resolution 
Picture Transmission (HRPT) station telemetries include typically 4000 to 5000 scan 
lines of data (depending on local topographical/anthropogenic obstructions), translating 
to 11-13 minutes of capture time from the nominal initial-capture time stamp on the 
filename. 

[00281 Tho lOOSOl The primary function of SeaWiFS is to detect and monitor ocean 
color. Specifically, it has been designed to retrieve the small component (.about.Syo f-Syo 
of total signal) of water leaving reflected radiance whose spectral signatttre provides 
usefiil information about ocean phytoplankton content (expressed in terms of chlorophyll- 
A retrievals), A varying along-track sensor tilt angle of +/ 19.8.d e gr ee . +/- 19.8° is used 
to avoid regions of sun glint (specular reflection of the solar disk upon the ocean surface) 
where chlorophyll-A retrievals cannot be performed. The information obtained from 
such retrievals is valuable to marine research and industry interests. Applications include 
guidance to the fishing industry, diver visibility information, surface currents, physical 
oceanography, general marine biology, and climate impacts on the marine ecosystem. 
While atmospheric constituents such as cloud, dust, and aerosol are considered as noise 
to ocean color algorithms, the current research regards these components as the signal of 
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interest. This description of the invention report details how SeaWiFS data are applied to 
enhance atmospheric and surface features/phenomena. 



100301 Th e f00321 The main advantage provided by SeaWiFS over optical spectrum 
instruments flown on the geostationary (e.g., the Geostationary Operational 
Environmental Satellite (GOES) constellation, Meteosat 5,7, Geostationary 
Meteorological Satellite (GMS). and the hidian National Satellite System (INSAT)) and 
other polar orbiter (e.g., National Oceanographic and Atmospheric Administration 
(NOAA), Defense Meteorological Satellite Program (DMSP), and many others) 
platforms for the current application is its multiple narrowband visible channels that 
provide a true color imaging capability. 

[003 1 ] Instrum e n t [0033 1 Instrument Limitations and Other Caveats 

[00321 The |00341 The main deficiencies of SeaWiFS in terms of 24-hour operational 
utility are its lack of infrared channels, its limited temporal resolution, and proprietary 
nature. Infrared data are useful for nighttime imaging, cloud particle size estimation, and 
airborne dust enhancement over land surfaces, among other uses. At mid to low 
latitudes, only one to three passes are available near local noon providing at b e st a 
"snapshof' noon — ^providing at best a "snapshot'^ of the synoptic-scale situation on any 
given day. Captures are increasingly frequent at the higher latitudes owing to the orbital 
geometry. Typical captures for the regional centers currently processed are on the order 
of 2 to 3 per day. Thus, SeaWiFS products are most useful when juxtaposed with 
geostationary data that offer the temporal information required to infer details of the 
atmospheric circulation. 



Sciences Corporation (OSC), the data are fonnally available to the general public on a 
research-only basis and are provided at cost. The National Aeronautics and Space 
Administration (NASA) has contracted with OSC for the research use of a 5-year data set 



r00291 Inotrum c nt Adx antages [00311 



Instrument Advantases 




?[0035 1 Because SeaWiFS was built, launched, and operated by Orbital 
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commencing in September 1997. Researchers who have registered as authorized 
SeaWiFS user (r o giGtration at http://seavvifG.gGfc.naGa.gov/SEAWIFS.html) are granted 
permission to obtain the data from the NASA SeaWiFS archive, approved NASA Space 
Act Agreement ground station, or authorized individual. As authorized receiving 
stations, the NRCs providing near-real-time LAC data for this research are covered under 
the existing agreement, with the only criterion being that NRL does not disseminate the 
LAC-derived products to the general public in near real-time. 

[003 4 1 Data Proc c GGing 100361 Data Processing 

[00351 S e veral f0037| Several options exist for processing the SeaWiFS telemetry, 
depending on the format of the data. The formats currently received at NRL are NASA 
Level-0 (LO; 10-bit Data Capture Facility frame- formatter files), with Terascan Data 
Format (TDF) (pre-calibrated and navigated) upon special request. Data from the NRCs 
are in NASA LO format, and require several preprocessing steps to calibrate and navigate 
the data as described below. For reasons described further along, the current approach to 
creating SeaWiFS products is based on NASA LO data and is independent of the 
Terascan software package. Several approaches have been iterated upon in arriving at the 
ciurent procedure including various blended combinations of software packag e s - 
sem epackages — some more appropriate for rapid deployment to remote operational 
centers and others geared toward in-house processing with all the conveniences of locally 
available third party software. While the purpose of this d e acription report is to 




summarize the current approach, alternative options and their associated data/software 
requirements will be touched upon in passing. 



the Linux platform, was selected as the scripting language for automated processing of 
SeaWiFS data. Perl draws from the strengths of many popular scientific programming 
languages such as C++ and Fortran while enhancing on the capabilities of other standard 
shell scripting languages. The Perl scripts written in support of the current SeaWiFS 





tf0038] Worth noting is that tlie Perl programming language, running on 
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processing perform the requisite secretarial operations of data retrieval, file and directory 
manipulation, execution and status monitoring of telemetry updating and data processing 
codes, dissemination of imagery products to the website, and general housekeeping. Perl 
remains the lone constant among all processing options discussed in this report. 

r00371 S e aDAS Processing Option 100391 SeaDAS Processing Option 

1003 8 1 Ono l0040l One of the available SeaWiFS processing options is the GSFC 
SeaWiFS Data Analysis System (SeaDAS) software package (Fu, G., K. S. Baith, and C. 
R. McClain, 109 8 : SeaDAS: th e SoaWiFS data analysis system, Proc. of 'l.oup.th Pacific 
Oooan Romoto Sonoing Conf.. Qingdao, China. Jul. 2 8 31, 199 8 . 73 79). (Fu et a/., 
f 1998]). which was designed specifically for viewing and basic manipulation of SeaWiFS 
data. It uses an embedded runtime Interactive Data Language (IDL) license as its 
operating paradigm. The SeaDAS software is designed to operate on NASA LO, Lev e l - 1 
A (LI A; Level-lA (LIA; raw spacecraft and instrument telemetry retained with 
convertedgee- geo- location. instrument telemetry, and selected spacecraft telemetry), and 
Level-IB (LIB; sensor calibration applied to the data). LIA and LIB data are stored in 
Hierarchical Data Format (HDF). The software includes a graphical user interface 
enabling straightforward viewing and manipulation of the data, as well as command-line 
functions useful for incorporating into scripts as batch jobs. The SeaDAS package was 
used early on in this research to examine the general utility of SeaWiFS data in terms of 
coverage and data quality before proceeding with customized application. While 
SeaDAS is not used in the current processing (in the interest of gaining a higher level of 
control over product display options), several functions were adopted from this package 
for conversion from NASA LO to LIB HDF files and data correction. SeaDAS remains 
an extremely useful stand-alone tool for SeaWiFS data visualization, and is available via 
the Int e rn e t at http://s e adas.gsfo.nasa.gov. Th e websit e includes v e rsions visualization. 
Versions of SeaDAS for several operating platforms are available andhas a limited level 
of helpful online support, support exists. 
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>t0041| The Currently Implemented Processing Scheme 



fGG40fAI0042L 



_A simple overview of the current processing flow is as follows: I) 



SeaWiFS LAC telemetry captured at an NRC site and converted immediately to NASA 
LO formatted file. II) Data sent to the satellite data storage facility at Naval 
Oceanographic Center (NAVO) at Stennis Space Center in Mississippi via file transfer 
protocol (ftp) push, III) Data are pushed via ftp to NRL Monterey from NAVO, IV) A 
periodic local cron job searches for receipt of new files, and if found: V) SeaWiFS 
processing script commences to calibrate, clean, apply atmospheric corrections, and 
generate products, VI) Products are uploaded to a password protected NPRNET website. 
Latency of near real-time products is generally on the order of 1 .5 hours fi:om time of 
ingest at the NRC site. The remainder of this section provides additional details on the 
steps outlined above. 

rOOl n NRC Stations [00431 NRC Stations 

[00 4 21 The f00441 The mid-latitude NRCs included in the current NRL processing are 
Rota [36.63N, 6.34E]. Bahrain: [26.16N, 50.39E], and Yokosuka [35.28N, 139.67E]. 
Each station captures SeaWiFS passes on or around its local noontime hours (1100 - 
440 0(1100—1400 Z for Rota, 0700 1000 0700—1000 Z for Bahrain, and Qi^ 
04^ 0130—0430 Z for Yokosuka). NRL receives the NASA LO data from NAVO as 
Gnu-compressed files typically 1 hour after the nominal capture time. File sizes vary 
according to pass coverage, with typical values ranging from 20 to 80 Megabytes. The 
55 Gigabyte hard drive currently used as a storage bay for raw data, algorithms, and 
products is sufficient for several months of online storage. All data are eventually 
archived to CDROM or tape storage in a hierarchical structure partitioned by station, 
year, month and day. 



[00 44 1 A [0046] A cron job launched every half hour initiates the processing of 

newly received SeaWiFS data. A loop over all files present in the ftp "r e c e ive 




[00451 



Data Calibration 
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dir e ctorv" "receive directory" checks for any new and unprocessed data. If a new file has 
arrived, it is uncompressed and copied to the local storage drive for subsequent 
processing to fciB LlB calibrated radiances (expressed in mW/cm.sup.2 sr 
TgHhia^r mW/cm'-sr-iam). Binary executables taken from the SeaDAS software package 
are used for this purpose, and require an updated SeaStar ephemeris for correct earth 
location (SeaWiFS requirement of 1 -pixel geolocation accuracy). These data are 
available at ftp.samoa.gsfc.naGa.gov, and are updated on the local system if the current 
elements file is found to be more than one day old. Any errors incurred during the 
uncompressing or calibration steps are reported to a log file and the script is exited 
gracefiilly. 

[00151 Sccto r [00471 Sector Definition and Noise Corrections 

[00 4 6] Upon f00481 Upon completion of tlBLlB processing, an DDL script for 
customization of the data is spawned. Latitude, longitude, sun/sensor geometry, and 
SeaWiFS channel variables are extracted fi-om the HDP data. A table of pre-defined 
sectors for each station is included and readily updated for new operational areas of 
interest or one-time case studies. Definition of these sectors is based on the specification 
of latitude/longitude boxes (lower left and upper right comers). Owing to the scanning 
geometry, regions near the edge of the swath suffer spatial resolution degradation. As a 
result, imagery for these pre-defined sectors may include only partial coverage at 
compromised quality. A minimum data check of 20% coverage over the sector in 
question is enforced prior to proceeding with the processing. To ensure that the optimal 
information content for every SeaWiFS pass is provided, a "float e r" "floatef ' sector is 
included which selects a box data away from the swath edges (poor resolution) and the 
horizons (avoiding numerous scan line drop-outs). The trade off for high image quality 
is a random coverage area from pass to pass. 



[00 4 71 Once [00491 Once the sector of interest has been defined (or loops over 
numerous sectors), an embedded loop over SeaWiFS channels is performed for 
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corrections of pixel noise and line dropouts. Noisy pixels appear in the imagery as a 
^^sal t"salt and pepper e ff e cf effect" owing to the presence of spurious very high or very 
low data values. Line dropouts appear in imagery as black scan lines devoid of data or 
filled with noise. These imperfections are present in the data due to a combination of 
possible hardware limitations/problems and clutter associated with nearby structures. 
Cleaning of the data is accomplished in two ways. Lines identified as dropouts are 
replaced by previous good scan lines (bad lines at the very beginning or end of a data 
stream are replaced by their nearest good neighbor). If multiple adjacent dropouts are 
present, the previously patched line is used once again. This correction loses its cosmetic 
effectiveness for large data gaps, but provides marked improvements in cases of spurious 
line dropouts. Noisy pixels, assumed to be spurious, are identified using a threshold 
based on the absolute values of pixel derivatives computed along a scan line. WTien 
flagged, they are replaced by an average of the two adjacent pixels on the same scan line. 
Noise pixels on the edges of lines are simply replaced by the value of their nearest 
neighbor. FiG rFig. 2 illustrates the result of corrections applied to a subsector of 
SeaWiFS data collected over the Alps. The overall effect of these bad line and pixel 
noise corrections is dramatic, and regarded as an essential step in the production of high- 
quality SeaWiFS imagery. Additional noise corrections (not yet implemented) may be 
required to correct pixels where one of the channels is missing data (but was not flagged 
by the noise threshold). This results in a magenta, yellow, or cyan-biased pixel values in 
the RGB products. 

fOO4S4 f0050| Atmospheric Corrections 

[00491 Th e lOOS II The effects of atmospheric molecular (sometimes referred to as 
"Ravleigh") "Ravleigh") scattering must be removed from the data before creation of true 
color products. The amount of Rayleigh scatter is proportional to ^ (wavelength) and 
hence is strongest for the shorter (e.g., blue, violet) wavelengths. If left uncorrected it 
will manifest in the imagery as a semi-transparent milky blue haze increasing in opacity 
toward the edges of the swath (as higher sensor angles peer through an optically thicker 
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atmospheric path), detracting from image contrast and clarity. To remove this 
undesirable component of the satellite signal, we use a radiative transfer model (Mill e r, S. 
D.. G. L. Stephens, C. K. Drummond, A. K. Hoidingor, and P. T. Partain, 2000: A 
multio e nisor diagnostic cloud prop e rty r e tri e val schem e , J Geophys. Rss., 105, No. D15, 
1Q95 1997) (Miller et al.. [20001) to simulate the clear sky Rayleigh scatter over a dark 
surface as a function of the solar zenith, satellite zenith, and solar-satellite relative 
azimuth angles. Pre-calculated look-up tables were created for each SeaWiFS channel 
and the correction was applied on a pixel-by-pixel basis, provided the sun/satellite 
geometry in the Level IB data. The primary caveat associated with this correction is the 
assumption of a cloud-free Rayleigh contribution that in some instances may result in an 
over-correction of pixels containing high clouds near the edges of the swath. Visual 
inspection of imagery has not indicated a significant distortion owing to this effect, 
however. Regions near the edge of llie SeaWiFS swath should in general be interpreted 
with cautious regard to the possible distortions owing to larger atmospheric path lengths, 
three-dimensional effects (parallax), and possible bow-tie effects. 

-100521 Scaling 



[00511 Onc e [0053l Once the data have been calibrated, cleaned of noise, and corrected 
for Rayleigh scalier, the spectral radiances (I.sub.. lambda. ) £/^ are converted to 
equivalent isotropic reflectances (R.Gub..lambda.)£ g^ according to the relationship;- 
IoFo,,(l) 

^.=-^^ m 



[00521 wher e .mu..oub.o where [i n is the cosine of the solar zenith angle and F^sv^o is the 
solar spectral flux at wavelength .lambda ..!. By definition, the isotropic reflectance 
represents the reflectance produced by an equivalent lambcrtian surface, and varies from 
0.0 to 1.0. Any spurious values greater than 1.0 or less than 0.0 are thresholded by these 
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upper and lower bounds. The reflectances are scaled by /ogrsubvio in an attempt to 
equalize the imagery and prevent the very bright cloud/snow pixels from dominating the 
information content of the products. These scaled reflectances are byte-scaled over the 
256-element color table range of the IDL palette, with a minimum logarithmic reflectance 
of 1.65 1.65 mapped to the minimum saturation value (0) and 0.0 mapped to the 
maximum (255). This proc e dur e is d e tail e d in th e comput e r program code lioting in the 
App e ndix. 

[00531 Tho [0054l The scaled data are warped to a mercator projection corresponding 
to the sector of interest. This procedure is currently the most computationally time- 
consuming component of the processing, and options for optimization or alternative 
approaches are currently being explored. 

[005 4 1 Tho lOOSSl The final image products are created as jpegs at 80% of full 
quality. File sizes range from 200 to 800 Kilobytes depending on the sector and image 
information content. The true color, vegetation enhancement, and dust enhancement 
products for various stations and subsectors are copied to a password-protected online 
holding directory on the NRL Satellite Meteorology web page 
(http://lcauai.nrlmry.navy.mil/sat producto.html). page. 

[00551 Tru e Color Enhancem e nt 100561 True Color Enhancement 

[00561 To [0057l To produce the true color images, the blue (412 nm), green (555 
nm), and red (670 nm) SeaWiFS data are loaded into the respective blue, green, and red 
IDL "color guns." "color guns." The data comprise a single matrix of dimension (n,m,3) 
where n=along-scan direction and m=along-track direction. The resultant imagery is 
similar, but not identical to what would be observed by the human eye from the vantage 
point of outer space; clouds as shades of gray to white, ocean surfaces as shades of blue, 
and land features as earth-tones and gi-een in regions of heavier vegetation. 

[00571 Examplcs [00581 Examples of true color SeaWiFS imagery are shown for the 
Rota, Bahrain, and Yokosuka regional centers in HGSr Figs. 3-5, respectively. FIG ^Fig. 3 
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demonstrates the immediate ease with which significant airborne dust is identified over 
oceans. The image depicts a widespread easterly plume of Saharan dust captured within 
the synoptic flow of a baroclinic system off the coast of Spain and Portugal. SeaWiFS is 
useful in monitoring these dust plumes as they migrate across the Atlantic Ocean (often 
reaching the eastern United States). The green colors in littoral regions of the Persian 
Gulf in HG rFig. 4 correspond to a higher amount of bioactivity (e.g., plankton) in 
nutrient-rich waters. Cloud streets over land and evidence of biomass burning in the 
form of diffuse smoke plumes are apparent over Southeastern Asia in FIG rFig. 5. The 
NRG in Yokosuka also captures seasonal dust storms blowing off the Gobi desert. Owing 
to the general circulation of the atmosphere, large Gobi and Saharan dust storms may 
reach the continental United States (the fomier crossing the Pacific with the "Roaring 
Fourti e s" ; Roaring Fourties". and the latter crossing the Atlantic with the Tropical 
Trad e s. Trades). 



[005 8 1 Vegetation Enhancomont [00591 Vegetation Enhancement 



[00591 Bccaus e [0060| Because green vegetation possesses a relatively high albedo 
compared to soil in the "refl e ctiv e " "reflective" or shortwave infrared (0.7-1.3 
■mu.m)v |im), a combined vis/shortwave-IR product that enhances vegetation regions is 
generated in a way similar to the true color enhancement described above. This 
vegetation enhancement product is similar to what has been produced for many years for 
land-usage studies with Landsat TM data. The 670 nm (red) channel in the previously 
discussed RGB true color enhancement is replaced by the 865 nm (SeaWiFS channel 8) 
channel data. This results in green vegetation appearing as bright red and appears usually 
in sharp contrast to surrounding regions devoid of vegetation. An example of this 
product is shown in FIGrFig. 6. The vegetation enhancement may prove useful in 
identifying potential fetch areas for airborne dust. 

[00601 Airborne Duot Enhancement 
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[0061] Airborne Dust Enhancement 

f006n An f 00621 An important capability of SeaWiFS, resulting from its multiple 
narrowband visible channels and indicated in the examples above, is discrimination 
betv^een airborne dust and clouds/smoke via tme color methods. Fine-grain particles 
lifted from sparsely vegetated and dry desert surfaces by strong frictional winds (typically 
greater than 10 m/s, usually by meso- and s>Tioptic-scale circulation features) give rise to 
significant plumes of atmospheric dust. These plumes follow the low to mid-level 
atmospheric circulation and may also become incorporated within synoptic scale 
baroclinic systems (giving rise on occasion to unusually "mud "mud i=aifi ^rain" storms that 
have been reported to deposit a coating upon everything in their path). 

[00621 P e rhaps f00631 Perhaps more importantly to Navy interests, airborne dust has 
significant detrimental impacts on slant-range visibilities (with visible optical depths that 
can exceed 3.0 and accompanied by visibilities less than 1 nautical mile) and also poses a 
serious hazard in terms of its potential to damage turbine engines. This is of particular 
relevance to aircraft carrier operations in regions of the world prone to widespread dust 
outbreaks. Because airborne dust retains the earth tone hues (usually light brown or 
yellow, depending on the source region) of the underlying surface, it is readily 
identifiable in most cases from a true color image, especially for optically thicker dust 
plumes. For regions of less dust (where its effects on slant range visibility may still be 
appreciable), an additional enhancement beyond true color is required. NRL Monterey 
has developed a dust enhancement algorithm for SeaWiFS that attempts to provide this 
capability. Enhancement of dust is predicated on basic spectral properties that 
characterize its appearance at visible wavelengths. As indicated in FiG rFig. 7, yellow 
light (characteristic of many desert dust plumes) results from the removal of blue light 
(via preferential absorption by the dust). This spectral contrast is even more dramatic 
between the blue and shortwave infrared channels. This preferential absorption does not 
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occur to the same magnitude in liquid/ice clouds (explaining why clouds appear as shades 
of gray to white, indicating nearly equal reflected intensities of red/green^lue light, in a 
true color composite). By defining a dust enhancement parameter ".DELTA.": 2 = 8 65 
112 8 65 4 12 "A": 



[00631 whore ■alpha..sub..lambdQ. where a> are reflectances at wavelength .lambdaA 
(nm), any dust component of the image will be enhanced. Cloud pixel values, which are 
relatively spectrally flat, will be diminished in brightness owing to the small numerator 
difference. The logarithm of the .DELTA. A parameter is scaled between 0. 4 5 0.45 and 
0.20 and loaded into the red channel of the RGB composite. 

[006 4 1 Ovvin g f00641 Owing to the significantly reduced red component for cloud pixels 
in the dust enhancement (small values of the .DELTA. A parameter) image, clouds 
lyi^ically appear as shades of light blue or cyan. Retaining the original grayAvhite cloud 
hues in the dust enhancement (for cosmetic purposes only) requires the implementation 
of a simple cloud mask - pixels mask — pixels flagged as cloudy can be re-assigned to their 
original red channel values. Clouds in the imagery are identified first by computing the 
mean scaled magnitude and standard deviation of the 412, 555, and 670 nm channels. 
Mean magnitudes exceeding 50% and standard deviation less than 2.5% are used as 
thresholds for positive cloud identification. The simple premise behind choosing these 
thresholds is that clouds are in general bright (higli brightness magnitude) and relatively 
spectrally flat over the visible (low standard deviation). Pixels flagged as 
"oloudv" "cloudv" are reset to the red channel value in the AA parameter. As the dust 
enhancement enhances land surfaces also, a land mask is applied to this product. 
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[00651 FIG. [00651 Fig. 8 demonstrates the dust enhancement product covering the 
Persian Gulf region. The top panel pair reveals that while the true color imagery is useful 
in identifying the significant dust plumes, the dust enhancement reveals a far more 
extensive coverage of airborne dust circulating within the Gulf. The center and lower 
panel pairs further illustrate this enhancement capability. While not necessarily posing 
the same hazard to aircraft engines as the optically thick (bright) plumes, these areas of 
enhanced aerosol have detrimental impacts on slant-range visibility. Although more 
difficult to detect than over the oceans, the presence of significant dust over bright and/or 
sparsely vegetated land surfaces can in some cases be inferred by noting contrast 
reduction in the background terrain. The lack of an infrared channel on the SeaWiFS 
instrument precludes the technique of detecting elevated (cooler) dust above a warm 
surface background via emission measurements. The utility of EOS-MODIS, which 
provides 36 channels (including narrowband red, green blue, NIR and LWIR bands), will 
be explored in this capacity. 

[00661 OthG r [0066| Other Utilities 

[00671 Th e [00671 The global coverage afforded by the polar orbiting platform 
provides SeaWiFS with opportunities to sur\'ey many atmospheric and terrestrial 
phenomena. Its utility as an additional sensor for tropical storm monitoring (e.g., 
Typhoon Pabuk showTi in FiGrFig. 9) may increase in relevance in light of the current 
mechanical problems aboard the Geostationary Meteorological Satellite (currently 
providing coverage of the Tropical Western Pacific and the primary satellite observation 
platform supporting the Joint Typhoon Warning Center). Unexpected events provide 
unique opportunities to examine additional capabilities of the SeaWiFS channels. The 
eruption of Mt. Etna, shown in HG rFig. 10, demonstrates an additional utility of the dust 
enhancement product r e vealin e product — revealing features of the ash plume difficuh or 
impossible to detect in panchromatic or even true color renditions. Depending on their 
composition and concentration (not all volcanic plumes are comprised of materials that 
will be enhanced by the currently detailed dust enhancement), such plumes may be 
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tracked around the globe, providing invaluable guidance to the commercial airline 
industry for avoiding these potentially serious flight hazards. The detection capabilities 
of SeaWiFS offer a valuable asset to Volcanic Ash Advisory Centers (VAACG; e .g., 
http://w>\w.oad.noaa.govA^AAC/washington.html) distributed world-wide whose charter 
is to monitor and track these plumes globally, and also the Geneva-based World 
Meteorological Organization (WMO; e.g., http:/A v ^v.wTno.ch/indoxflash.html) , which 
coordinates global scientific activity among an international panel of participants. 

[ 006 8 ] - Now Proc e ssing Formats [00681 New Processing Formats 

[00691 Th e i00691 The SeaWiFS processing as discussed in this description r eport is 
well suited for the creation and analysis of high quality SeaWiFS imagery products in an 
automated, near real-time framework. 

[0070] For l0070| For the reasons cited above, efforts to reproduce this processing in 
a fashion entirely independent of third party software have resulted in an analogous 
package that assumes receipt of SeaWiFS data in TeraScan Data Format and proceeds to 
customize the data in the FORTRAN-90 programming language. True color and dust 
enhancement image rendering is then completed using the "rgbimag e " "rgbiniage'^ 
TeraScan intrinsic function. This FORTRAN-based version is actually more efficient 
and hence runs faster than its IDL counterpart. 

[00711 Blondod Producta [0071 1 Blended Products 

[0072] S e aWiFS f00721 SeaWiFS is only one of several telemetries available for 

processing at NRL Monterey. It is a straightforward procedure to incorporate data from 
these other sensors within the same image to produce a multi-sensor product. For 
example, the masked-out land areas in the dust enhancement product may be filled in by 
temporally matched geostationary data (e.g., Meteosat 5 over Bahrain, MeteoSat 7 over 
Rota, and CMS over Yokosuka). Infi-ared information can be used to enhance regions of 
elevated dust over the land in certain situations. As NRL also has at its disposal analysis 
and forecast data from its mesoscale (Coupled Ocean-Atmosphere Mesoscale Prediction 



Inventor: Steven D. Miller 



Serial Number: 10/713.908 



PATENT APPLICATION 



Navy Case 84.024 



Svstem').TM. Svstem)^" and global scale (Navy Operation Global Atmospheric Prediction 
System) models, it may also prove useful to blend these data into the products. Using 
again the example of the dust enhancement product, inclusion of the 
COAMPS.TM. COAMPS^" 925 mb wind field in the form of vectors overiaid upon the 
imagery may serve not only as a qualitative validation of the model analysis but also 
provide the METOC officer in the field with a means to predicting the migratory 
behavior of significant dust plumes in the short term. Overlay of short-tenn forecasts of 
the frictional wind over land (winds near the surface that are responsible for lifting the 
dust into the atmosphere) may also be of use. Of course, any such multi-sensor and 
model-fusion products will require a training module identifying the various components, 
strengths, and weaknesses. 



validation of the Navy Aerosol Analysis and Prediction System (NAAPS; a global, multi- 
component analysis and modeling capability to produce aerosol products from satellite 
data and surface-based measurements). The goal of NAAPS is to improve the 
forecasting/nowcasting of visibility restrictions caused by aerosols, including dust, 
smoke, and volcanic ash. Information and for e cast produota associated with NAAPS may 
bo found at http://www.nrbiiry.navy.mil/aerosol. 

[0071] Extension to MODIS (Moderate Ro o olution Imaging Spoctroradiomotor 
Proceaning') [00741 Extension to MODIS (moderate resolution imagine 
spectroradiometer) Processing 

f00751 Tho r0075l The MODIS instrument offers a similar true color imaging 
capability together with improved spatial resolution (up to a factor of 4 over SeaWiFS 
HRPT). HGr Fig. 1 1 provides an example of the MODIS dust enhancement over the 
Northern Arabian Sea. Plumes of dust over ocean and dry lakebeds on the desert interior 
(often serving as sources for dust) are enlianced as shades of pink. Of specific relevance 
to the dust enhancement product is the availability of thermal infrared channels that 




data are of immediate relevance to both the initialization and 
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enable the discrimination of elevated (cooler) dust plumes over warm backgrounds. See 
the Appendix for th e computer program code listing implomonting thist e chniqu e . On- 
board calibration of MODIS radiances supports accurate physical retrievals of aerosol 
optical properties. A technique for producing 250-m true color imagery using 500-m blue/green 
and 250-m red channels has also been developed. This high spatial resolution capability is 
enhanced by temporal coverage when considering the Terra (MODIS: 1030 descending node). 
SeaStar (SeaWiFS: 1200 descending node), and Aqua (MODIS: 1330 ascending node: platform 
to be launched in April 2002) platforms in concert. This quasi-geostationary 3-hour loop will 
provide valuable information on cloud and/or dust motion. An enabling docum e nt motion. 

for th e MODIS lond/ocoan sigDificant dust e nhancement will b e publish e d in the 
G e ophysical Research L e tt e rs. 



[00761 A [0077l A new processing capability for generation of true color, 
vegetation enhancement, and ov e r water dust enhancement imagery from the 8-channel 

from Navy Regional C e nt e rs in Rota, Bahrain, and Yokosuka currently populate the 
following password protoctcd w'ebsit e on developed by the NRL Monterey sat e llit e pag e 



A pplications Section team. 

[0077] Tw^o |0G781 Two main points become clear from the above commentary: 1) 
SeaWiFS and other sensors offering true color imaging capability are of significant value 
and in high demand for Navy operations, and 2) in spite of the fact that such data has 
been available for some time now from the satellite platform (e.g., LandSat-TM) this 
capability has not been transitioned effectively to the end-users who require it (i.e.. those 
deployed in the field). In turn, the twofold intent of this description report has been to 
demonstrate the current in-house processing capabilities for the SeaWiFS sensor and 
articulate this operational demand. A logical extension to these statements is that there 




-100761 



Conclusion 
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exists a pressing need for X-band receiving stations in perennially sensitive regions of the 
world for the purpose of capturing MODIS and other high-data-rate (both multi- and 
hyper-spectral) satellite telemetries. The potential savings in terms of both loss of life 
and equipment (R. Contu; "Tli e Rol e of W e ather in Major Naval Aviation Mishaps", MS 
Thosis. 106 pp.. NPS Monterey) which is rcportGd at (e.g.. Cantu 120011 reports a 51 
million dollars annual loss in Navy equipment attributed to visibility-related 
meteorological phenomena such as aerosols, dust, and atmospheric trapping layers) 
owing to improved METOC guidance gained from these state-of-the-art sensors far 
outweighs the one-time cost of the receiving stations themselves. 

[007 8 ] Although this inv e ntion has b e en describ e d in relation to on ex e mplary 
embodim e nt ther e of, it will b e und e rstood by thos e skill e d in the art that still other 
variations and modifications can b e aff e ct e d in the pr e f e rred embodiment without 
detracting from the scope and spirit of the inv e ntion as described in th e claims. 



